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ABSTRACT 


This sensitivity analysis employs mathematical models which estimate the total 
dose. Maximum ERD, and probability of casualty or fatality for an individual exposed 
to a particular radiological environment. The objective of the analysis is to 
determine the relative importance of the parameters normally employed in evaluating 
Civil Defense systems; i.e., to determine those parameters which, through large 
variance or inaccurate estimates, will contribute most to erroneous evaluations of 
CD systems. The sensitivity of the dose (or probability of casualty) to variations 
in the input parameters defining the radiological environment is examined. The 
fotal dose model is analytical and the ERD model is computerized. Ihe parameters 
examined are reference intensity, time of arrival, time outside in fallout, radia¬ 
tion decay exponent, ERD recovery fraction, ERD recovery rate, duration of fallout 
buildup, and protection factor. Sensitivity indices are calculated for each parame¬ 
ter. The sensitivity index is defined as the fractional change in dose (or proba¬ 
bility of casualty or fatality) divided by the corresponding fractional change in 
the input variable. It is concluded that dose and casualty computations are quite 
sensitive to errors in the field decay exponent, and they remain sensitive over 
the examined range. Sensitivity to variations in fallout reference intensity 
and protection factor are high over the whole range of parameter values. Sensitivity 
of time of arrival of fallout can be quite high for the lower values of the parame¬ 
ter. Sensitivity of dose and casualty computations to the remaining parameters is 
low in most cases of interest. Expansion of the sensitivity analysis to include 
parameters other than fallout, which define the total casualties from a given attack 
on the United States, is necessary before further conclusions concerning a national 
vulnerability analysis can be drawn. 
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A Gonoralizcd Scnf-itivity Analysis of CD ■Systcrr’^ 


I. INTRODUCTION 

A. Purpose 

An analysis of the effectiveness of the civil defense shelter system nust take 
into account many variables. These variables include some that determine the fall¬ 
out environment (fallout reference intensity, its time of arrival, its rate of build¬ 
up, and the radioactivity decay rate); some that determine the protection against 
radioactivity (the gamma ray shielding characteristics of available shelter, and the 
time required for travel through fallout before shelter is reached); and sor.e that 
Influence physiological condition (recovery rate, and the fraction of radiation 
damage that can be repaired). The research reported in this volume is directed at 
evaluating the importance of uncertainty or inaccuracy in each of these variables 
in estimating effectiveness of a shelter system in a fallout environment. Effective¬ 
ness is measured in terms of Equivalent Residual Dose (ERD), probability of fatality, 
and time required in primary shelter for the ERD to return to a specified level. 

The results from such an analysis can be applied in two ways: (1) determination 
of research priorities and requirements for increasing the validity of shelter 
systems analysis, and (2) improvement of the models or formulations that arc used for 
systems analysis. 

B. Procedure 

."1. General " - --- 

The procedure followed in this volume is known as sensitivity analysis, 
which is a systematic determination of the effects of small changes in input 
variables. Such an analysis requires the utilization of an equation or a 
mathematical model that relates the input and output varlables--such as reference 
intensity and estimates of casualties. Each input variable is varied, whilt 
all others are kept constant, so that the effects of these variations on the 
output are measured and compared. 
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// 


Specific 



Because of the number and nature of the variables that Influence estimates 
of the effects of the radioactive fallout, Che establishment of a model for 
sensitivity analysis Is cor^lex. Ihe major steps that were followed in this 
analysis are given below and are covered in more detail In Section II, which 
follows. A generalized flow chart of the procedure Is shown in Figure 1, 
Research Steps in Sensitivity Analysis. 

Step I. Selection of Parameters and Values . 

Several values for each of eight factors contributing to radiological 
effects were selected as input parameters. Combinations of these values were 
used later in equations (derived while performing Step II) to provide the cases 
upon which the analysis was based. (See Section II.A below.) 

Step II. Derivation of Basic Equations to Represent Fallout Effects . 


Equations were derived for determining ERD and corresponding probability 
of fatality (Pj^) on the basis of the eight selected factors. (See Section 
II.B below.) 

Step III. Selection of Formulas for Determining Relative Sensitivity of 
Input Variables . 

Ihe following relationship was established as the sensitivity index: 





Where S^ is the sensitivity index; Y is the measure of effectiveness for a ' 
selected set of cases; aY is an arbitrary value change in Y; is the value 
of any single input parameter in the same set of cases; and aX^ is the change 
in this input necessary to cause aY. (See Section II.C below.) 


Although it would have been desirable to relate the equations in Step II 
directly to the sensitivity ratio in Step III, it was not possible to do this 
readily without making over-simplifying assumptions. Consequently, intermediate 
steps were required. 


Step IV. Compilation of Cases for Analysis. 

Several levels of each input variable were selected’ and the Maximum ERD 
and probability of fatality were computed for each combination of input 
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variables. Using a CDC 1604 computer and the Mainline Program which was designed 
specifically for this purpose, approximately 43,000 cases were developed. After 
those cases were excluded in which the Maximum ERD was less than lOOr or greater 
than 700r, a total of 10,143 cases were left for the sensitivity analysis. 

(See Section II.D below.) / 

Step V. Grouping the Cases and Calculation of Mean Values . j 

' A tape containing the 10,143 cases was transferred to a general purpose 
data processing program written for the IBM 7072 [Reference 1] in order to group 
the selected cases into subsets for analysis. (See Section II.E below.) 

Step VI. Determination of Model Approximation Equation . 

Using the 10,143 selected cases and multiple regression, the same general 
purpose program [Reference 1] was utilized to determine a single equation 
(the principal model approximation) fur the calculation of Maximum ERD as a 
function of the eight input variables. (See Section II.F below.) 

Step VII. Determination of Principal Sensitivity Indices . 

The effect of each input parameter on estimating fatalities was calculated 
using the sensitivity formula (Equation 1) stated in Step III, the mean values 
of the input variables determined in Step V, and the principal model approxima¬ 
tion referred to in Step VI. (See Section II.G below.) 

Step VIII. Verification of the Principal Sensitivity Indices . 

The change in sensitivity indices was examined throughout a range of values 
using a Total Dose Model. (See Section II.H below.) These results were com¬ 
pared to those contained using the model approximation equation of Step VI. 

(See Section il.I below.) 

Step IX. Sensitivity Analysis of the Effect of Input Parameters on Shelter 
Stay Time . 

A subgroup of cases was used in an analysis of sensitivity in which the 
dependent variable was the shelter stay time required for the ERD to return 
to 80r. (See Section II.J below.) 

Step X. Evaluation of Significance of Individual Sensitivity Indices . 

The sensitivity indices of the various parameters were then evaluated in 
view of the accuracy with which these parameters could be expected to be 
established as inputs to any analysis of system effectiveness. (See Section 
II.K below.) 
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II. THE SENSITIVITY ANALYSIS 


A. Selection of Paranctcrs and Values 

In order to obtain a representative sample of cases on which to base the analysis 
of sensitivity, a range of reasonable values was assigned to each of the input variable 
commonly used in shelter systems analysis. These values are summarized in Table I. It 
will be noted that in all cases except that for protection factor (PF), Che values are 
evenly spaced. 


TABLE I 

Values and Symbols for the Input Parameters 
Used in the Sensitivity Analysis 


Parameter 

Symbol 

Values 

Reference Intensity (r/hr) 

I 

o 

300, 1500, 2700, 3900 

Time of Arrival (hrs) 


1, 4, 7 

Buildup Duration Factor 

E 

1.13, 3.64, 6.15 

Protection Factor 

PF 

2, 10, 40, 100, 500 

Time in Fallout (hrs) 

^2 

0.0, 0.3, 0.6 

Field Decay Exponent 

Z 

-1.0, -1.2, -1.4 

ERD Recovery Fraction 

F 

0.85, 0.90, 0.95 

ERD Recovery Rate (Fraction/day) 

B 

0.020, 0.025, 0.030 


Fallout hazard (Risk) data on 200 Standard Metropolitan Statistical Areas 
(SMSA's) were used as the basis for selection of values of reference intensity and 
time of arrival. These data were extracted from an OCD source. Times of arrival 
between one hour and seven hours occurred in about 50 percent of those cases; 
reference intensities were between 300 r/hr and 3900 r/hr for 94 percent of the cases. 
(See Volume II, Figs. 31 and 32.) 

The buildup duration factor (E) is a factor by which time of ariival is multiplied 
in order to represent estimates of possible time between the beginning and cessation 
of fallout. Appendix A contains a discussion of buildup functions. 

The protection factor (PF) values were based partially on National Fallout Shelter 
Survey data summarized in Reference [3]. Identified shelter spaces with PF's of 40 
and above could accommodate roughly 20 percent of the United States population. PF's 
of 40, 100, and 500 were selected to represent the three levels of PF's between 40 
and 99, between 100 and 250, and greater than 250, respectively. Factors of 10 ind 2 
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were used as representatives of the PF of population in residential basements and in 
homes without basements, respectively, , 

"ime in fallout, prior to shelter occupancy was set arbitrarily at 0.0, 0.3, 
and 0.6 hours. 

In determining equivalent residual dose (ERD), values used for the radiation 
field decay exponent and for the ERD recovery fraction and ERD recovery rate were 
varied about the generally accepted values used in systems analysis. 


B. Model Basic Equations 

The following basic equations constitute the analytical formulation of the 
“mbdel.^ It was upon these equations that the Mainline Computer Program (see Appendix 
B) was developed.. 

1. Equivalent Residual Dose (ERD) 


ERD (n) - D + d 
' ' n n 


I>n - a-B) Vl + ^ /n-l 

‘'n - ‘^n-l 


( 1 ) 

( 2 ) 

( 3 ) 


where: 


D ■ Recoverable dose remaining during time interval n-1 to n 
n 


dose remaining during time interval n-1 to n 

I 

,1 

a of recovery for reparable damage (fraction/ 


(roentgens). 

d ■ NonrecoverablJ 
n 

(roentgens). 

B B Fractional rat] 

--unit time). 

I 

F >■ Fraction of dose which is reparable. 

(1-F) « Fraction of dose which is nonreparable. 

n « Number of time units after detonation, n * 1, 2 , 3, ..., 
when Increments (n) in time are 24-hour days. (Time is 

incremented in minutes for the first 24 hours and in days 

thereafter.) 

t ■ Number of hours after detonation. 

$(t) ■ Field Intensity in r/hr at time (t) after detonation. 


and: 


t(t) - 0 


for: t < T, 
““ A 










where: 











sensitivity index for each input parameter. 


We begin the mathematical development by first defining the sensitivity index 
as the ratio of the percent change in the output, or dependent variable, to the per 
cent change in an input parameter [Reference 4]. 

Thus, the sensitivity Index (S^) of the i^^ input parameter, measured from 
mean values, is: 

[Q. 'ge in the dependent variable (AY)] 

_ (Kean / ’ue in the dependent variable (Y)] 

. — . . . 

[Change in the input parameter (AX^^)] 

[Mean value of the i^^ input parameter 


or; 



( 4 ) 


D, Calculation of ERD and Probability of Fatality 

The model basic equations were used with each combination of input values in 
the Mainline Program^^ on the CDC 1604 computer. 


For each combination of input parameters, the time required for the buildup 
of the radiation field is calculated. Next, the dose a person would receive, if 
any, while outside the shelter and the dose at the end of the first day Inside is 
computed and totaled. Ihe dose is then split into a recoverable S'^d nonrecoverable 
fraction and a test is made to determine if the probability of fatality equals 

"one.” If so, the program returns to the initial state and selects another inpu^_ 

combi .tion. 

Itie ERD is accumulated and updated by recoverable and nonrecoverable fractions 
and the following tests are made on the daily increments: 

(1) If a Maximum ERD has been reached, its values and date of occurrence are 
recorded for output. 


— An additional independent variable, D,, was included to set three possible levels 
to which the ERD could return—80r, 50r, and 20r. There was a total of 42,740 
cases for all combinations of input values. 
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/ 

/ 

(2) If a Maximum ERD has been reached and the nonrecoverable ERD has not ex« 
ceeded the allowable dose, and the accumulated ERD Is less than the 
allowable ERD, then the day Is recorded for output. 

Eventually the test conditions will be met and the program will branch out of 
the day-by-day calculation of the ERD. Once out of this routine, the program will 
take the Maximum ERD value and calculate the probability of fatality and other de¬ 
pendent variables. ___ 

The dependent variables calculated by this program are; 

(1) Maximum ERD. 

(2) Probability of fatality. 

(3) Probability of being a noneffective. 

(4) Day in which the Maximum ERD occurs. 

(5) Day in which allowable ERD occurs. 

(6) Nonrecoverable dose on that day. 

Cases were excluded in which the Maximum ERD was grearer than 700r cr less than 
lOOr, thereby leaving for analysis a group of 10,143 cases (Gioup 0 of Table II, 
Basic and Subgroup Divisions of the Mainline Program Cases) . Thus, the sensitivity 
analysis was based on cases in which the interaction of input parameter values re¬ 
sulted ifl situations of maximum interest--that is, cases in which even small varia¬ 
tions in input parameters could affect both casualties and fatalities. 

E. Grouping of Cases and Calculations of Mean Values 

Subgroups of Group 0 (which contained all cases between lOOr anc 700r Maximum 
ERD) were established for each of the four input values of reference intensity. In 
addition, it was desirable to combine two of these subgroups in order to have a 
subgroup with "moderately” high reference intensity (Subgroup 20 of Table II). 

In order to examine the significance of the time-in-shelter data, another 
basic group of data was established (Group 1). Into this group were put cases in 
which the ERD would eventually reach 80r or less, because 80r is considered to be 
the threshold below which physical disability is unlikely to exist. Within this 
latter basic group, subgroups based on reference intensity, allowable ERD, and 
protection factor were established. Table II shows the basic and subgroup divisions. 

The cases were transferred from the Mainline Program by a tape to the IBM 7072 
computer. A general purpose data processing program (TSAR) at Duke University was 
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TABLE II 


Baste and Subgroup Divisions of the Mainline Program Cases 


BASIC GROUPS SUB- 

■GROUPS 

NUMBER 

SORTING PARAMETER 




OF CASES 


Group 0 



10,143 

Reference Intensity (r/hr) 

(parameter Combinations 
with : 

5 


2,115 

1-300 

lOOr < Max. ERD <" 700r) 

6 


3,108 

1° - 1500 


7 


3,431 

1° - 2700 


8 


1,48? 

1° - 3900 


20 

of 

(sum 

7 & 8) 

4,920 

1° - 3063 
o 

Group 1 



3,481 


(Cases which return to 
an allowable ERD < 80r) 




Reference Intensity (r/hr) 


9 


909 

I - 300 


10 


1,041 

1° - 1500 


11 


1,058 

1° - 2700 


12 


473 

1° » 3900 

0 





Allowable ERD (r) 


2 


2,069 

D - 80 


3 


1,204 

dJ - 50 


4 


208 

dJ - 20 





Protection 


13 


987 

PF - 2 


14 


1,291 

PF - 10 


15 


795 

PF - 40 


16 


321 

PF - 100 

■ 

17 


87 

PF - 500 





Multiple Sorting Parameters 


18 


352 

D^ - 80r 





0 < P- < .5 





2 — 





T 2 - 0 


19 


240 

D, - 50t 

4 





0 < P- < .5 





2 — 





T 2 - 0 



- 10 









used for the grouping of data into srts 'Rof 


The same program was used for determining the mean values of the input parameters 


for the basic groups of cases and for each of] 
of the Basic Parameters , shows the pertinent 
groups. 


Determination of the Principal Itodel Appiroxination 


fence 1 


the subgroups. Table III, Mean V^alues 
hcans for each of these groups and sub- 


For purposes of computation it was desirable to have a single equation that 

del. Using a stepwise regression to 
, logarithmic, and exponential) were 
fitted to the sets of data produced by that miodel. This was accomplished also by 


approximates the Mainline Program computer rro 
obtain this, three basic functions (quadratic 


the TSAR general purpose program. Each funct 


«[Reference 1’. Iflquation 5 best fitted the dajta and was made the principal model 
approximation 


lion was fitted by th r step procedure 


f(X) » Max. ERD “ 

^(-2.7468) j (.932 + .0001499PF 


- .0446F + 5.354B - .006493E) 


_ (.4346 - .7954Z) p_(-1.0949 + .041811, ¥ .06268Z) 
_(-3.074 - .1161E) „(-1.1553) j 


( 5 ) 


G. 


I 

The statistical characteristics of the principal model approximation equation 


are; 


(1) Multiple correlation coefficient; 


10.9991. 


(2) Stand'ard error of the estimate: 0.'02421 (i.e., c of ERD + 6r). 


(3) Value of F-ratio: 886.785 (a further indication of the high significance 
of the equation). 

(4) Data Group 0 (Table II). 


Calculation of the Sensitivity Indices i 

We now have a formula (Equation 5) relatling the input variables to the dependent 
variable Maximum ERD, and a formula for the sjensitivity index (Equation 4). To 
calculate the sensitivity indices at the mearjs. Equation 4 is placed in differential 
form; 


—^ Identification of the variables is in Tabte I. 
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Individuals in 500 PF shelters only reached or exceeded the MaxiTnum ERD if they receive appreciable dose outside. 
Their time-outside average was twice any other group. | 









Let; 

so; 



Y - f(X) 


=1 


(^) 


( 6 ) 


Note that f(X) is the principal model approximation (See Equation 5); 

^2 ^3 ^4 ^5 

£(X) - K I T. PF Z B 

O A 


(7) 


where; 

. Cj .932 + .0001499PF - .0446F + 5.354B - .006493E 

- .4346 - .7954Z 

- -1.0949 + . 04 I 8 IT 2 + .06268Z 

- -3.074 - .1161E 

» -1.1553 
-2.7468 

K ■ e 


The partial derivatives of Equation 5 with respect to each independent variable 
(X^) are; 


C--1 C, C, C, C, 

» K C, I ^ T. ^ PF ^ Z ^ B ^ 
1 o A 

o 

^ffyy *'2”^ ^3 ^4 % 

ai^-KCjI^ PF^z‘b 

- K K I T pr'’"' 
aPF *^2 o ^ 

Pf‘^3 /4-‘ /5 

az 4 o A 


aim 

SF 


» K K 5 
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F - .9 
B - .025 
E - 2.57 


and using the Equation 4, we calculate the Sensitivity Indices (S^). 

TABLE IV 

Principal Sensitivity Indices 


Index 

Variable Denoted by Index Subscript 

Value 


Field Decay Exponent 

-4.22 

^PF 

Protection Factor 

-0.96 

^I 

o 

Reference Intensity 

+1,02 


Time of Arrival 

-0.50 

^F 

ERD Recovery Fraction 

-0.30 

^B 

ERD Recovery Rate 

-0.14 

^E 

Buildup Duration Factor 

-0.18 

^T 

^2 

Time in Fallout 

40.05 


It is upon these values that the conclusions and recommendations will be 
developed. 


The above sensitivity indices are obtained for Maximum ERO as the dependent 
variable. However, over the dose range where probability of casualty is neither 
0 nor 1, the relation between probability sensitivity (PS^) and dose sensitivity 
(S^) is given by the relation; 

a ' (9) 

where p probability of casualty = a • Dose + b. Hence, the relative rank of the 
S^ remains the same when evaluated at a given set of input variables. 

H. Total Dose Model 

The major computational difficulties in determining the above sensitivity 
indices arise because of the use of equivalent residual dose instead of total dose 
as a measure of biological effect. This section will describe a simpler model which 
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// 

/ 

f 

allows the calculation of the sensitivity indices analytically. In this model the 
dose rate buildup function is assumed to be linear, and total dose is used as the 
dependent parameter. 


A new parameter, T^, is Introduced in this approach. is the time period from 


entrance into shelter until the cutoff of the total dose calculations. Thus, 

(T^ + T 2 ]+ T^) represents the complete time period from time of detonation until com¬ 


pletion of the total dose calculation. Using the parameters previously defined (see 


/ 


Table I) and the new parameter T^, the total dose D is given by the expression: 


^ ^o 

®-pF 


(1+E)‘'E 


( 


(PF-l) 




<1 




(T^ + T 2 + (1+E) 




The sensitivity indices can then be calculated directly 

(PF-l) T 2 


, ^ "2 ^o "2 

s_ ■ m — ^ "j:!” 


■Ij - SI, 


D PF 


(1+E)‘E 


Z + (T^ + T2 + T^) 


s - . .lA 

% ^^A ° 


D^PF T 


Z-1 




2 (1+E) 


^ [ffF-l,««,(!) . (Z-l) E^ 


c ^ zi. *0 
PF ■ SPF ‘ D 2DPF 


2 ., 2 „ 2 . 


V <i+^) E 


(T, - E‘T/) 


(1-Z) [ 


+ T 2 + (1+E) 


I-E] 


0 0 


-z 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 
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(15) 



■ ?^E * D 


I E 
o 


D PF (1+E) 


Z+1 


n-PF)'l+E (1-^2) 2 

7J.1 -) * -) 


2T E^ 
A 


+ UdLfiZll ^i-z . 


(16) 


Z 

" -Z * D “ D PF 


t 


2 T 

T, ET 

(PF-1) ^ + In :T^(1+E)‘. iT^d+E); ^[-1] 

A J 




(T^ + T 2 + - [T^(l+E)] 


“1 


(1-Z) 

(ife |<^A ^ ^2 + ^ ^2 + ^o) 


- ;t^( 1+E);^‘^ In (T^(l+E) 


(17) 


The total dose sensitivity indices (S^^) can be readily converted to probability 
of casualty sensitivity indices (PS^), if desired. In the region of interest 
(probability neither zero nor one), the dose is related to probability of casualty 
(p) by a linear relation: 


p = aD + b. 


(18) 


Hence; 



X. 

1 

P 


^ 


(19) 


Only dose sensitivity will be discussed in the remainder of this analysis, since the 
relative rank o' both S and PS are the same for a given set of values of the inde¬ 
pendent variables. 


In Table V, Comparison of Sensitivity Indices , the values of the S^ obtained 
from Equations 7 through 13 by evaluation at the mean values of the variables listed 
in section G above are compared to the values of S^ listed in Table IV (calculated 
from the ERD model). 
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Index Variable Denoted by Subscript ERD Model Value 


'otal Dose Model Value 


Field Decay Exponent 


Protection Factor 


Reference Intensity 
Time of Arrival 

I 

Buildup Duration Factor 


-4.22 

-0.96 


+1.02 


-0.50 


-4.36 


-0.99 

+1.00 

-0.50 

-0.25 


Time in Fallout 


40.05 


40.02 


In Figures 2 thru 6 , the values of the sensitivity Indices are plotted as a 

function of the indexed variable. While the. rest of the independent variables are 

held constant at the mean values given previously, two values for FF (39.1 and 500) 

and two values of T^ (96 and 240 hours) are plotted for each jeurve. Figure 7 shows 

the increase in the sensitivity index for 7 ^ fallout) as Z (the decay rate) 

increases. Figure 8 shows the dose sensitivity to T as a function of T for the 

O j o 

two protection factors. This curve indicates that total dose is not particularly 

' 0 • ! 

sensitive to T^ when T^ is measured after 100-200 hours, espedially for high PF 
shelters. In fact, where Maximum ERD occurs within the first jfew days, it can be 
expected that the results from the total dose sensitivity analysis model will be 

essentially equivalent to the results from the ERD model. 

V ' , • ■ ■ 

I. Validity of Ranking the Sensitivity Indices 


In Tables IV and V, the sensitivity indices were evaluated for one set of input 

i 

variables. However, the relative values of these might change when evaluated for 
radiological environments other than the one defined by the mean values in section G. 
To investigate this problem, the following conclusions were drawn based on Figures 
2 through 8 and other calculations of sensitivity indices: 


(1) Sj , sensitivity index for reference intensity (I^)* M 1 over the entire 
range of values; thus, I^ remains a relatively important parameter in all 


(2) S„_ is essentially constant at 1 for most cases of Interest; thus, protec- 
PF 

tlon factor (PF) is a relatively important parameter in all cases. (See 
Figure 2). 
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Protection Factor (pp) 

Fig. 2. Sensitivity Index of PF as a Functian of PF 
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Fig. 4. Sensitivity Index (S_) as a Function of Buildup Factor (E) 





































































Sensitivity Index (S^ ) as a Function of Time In Fallout 























































Fig. 6. Senaltlvity Index (S,) as a Function of Decay Constant (Z) 










Sensitivitj^JLndex for Time in Fallout (S ) as a Funct ion of Decay Constant 
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(3) S remains high for all cases, thus Z is extremely critical in vulnerability 
c$ 

analyses. (See Figure 6). 

(A) For analysis observing only Maximum ERD reached while in primary shelter, 

ERD concepts and a total dose approximation of Maximum ERD yield essentially 
the same results in sensitivity analysis. (This does not imply that they 
yield the same results in casualty calculations.) 

(5) Sj can become quite large for short times of arrival of fallout (T^) and 
hi^ protection factors. In areas where this combination of environments 

:_exists, must be treated carefully in vulnerability analyses. (See 

Figure 3). 

(6) In most cases of Interest, the other sensitivity indices remain appreciably 
less than 1, though some change significantly as the parameter values are 
varied. 

Analysis of Time Required in Shelter 
1. Procedure 

In the preceding analysis, the independent variable was Maximum ERD. It 
is also Important to know how the input variables may affect the time required 
in shelter after the beginning of fallout. 

In the analysis of time required in shelter, only those cases from the 
Mainline Program were considered in which the probability of a noneffective 
(P 2 ) was greater than zero and less than 0.5. Thus, the largest Maximum ERD 
to be considered was 200r and the smallest Maximum ERD was lOOr. (See the dose 
response relationships in Section II.B.2.) Another condition was that the ERD 
returned to 80r. It was assumed that the reduction of ERD to this level would 
mark approximately the end of the continuous stay in shelters. Time in fallout 
outside shelter was also excluded from consideration. 

The TSAR general purpose program was used to set up subgroups 18 and 19 
(see Table II) and to develop an equation by multiple regression. The equation 
given below determines the number of days in shelter until ERD returns to 80r. 

When: lOOr £ Maximum ERD < 200r, 








( 20 ) 


f(X)* - N, - 4.47 - 44Z + 2.26 ^ + 1.6D- - 1.55(Z x 

c r*. J rr 

+ 926(Z X B) - .7920(Dj x F) - 19.8(0^ x B) 

Multiple Correlation Coefficient: 0.9717 
Value of F-ratio: 830 


Where: 

N 2 ■ Duration of shelter phase until ERD falls to 80r 
Z ■ Field decay exponent 

^o « Reference intcnsity/protection factor 
PF 

“ Maximum ERD prior to N 2 

B ■ ERD recovery rate 
F ■ ERD recovery fraction 

The mean value of stay time for the cases considered (Y in the sensitivity 
formula) was found to be 47 days. 

The sensitivity indices were calculated using Equation 20, the differentia, 
form of Equation 4: 


a. 



c . ^f(X) 

i * 47 

witn; 

* Y * f(X) (see Equation 20), 

1^2 * Y = 47 days. 


( 21 ) 


2. Results 

The sensitivity indices for the input parameters are shown in Table VI, 
Sensitivity Analysis of Time Required in Shelter . It is seen that changes in 
the decay exponent (Z) are the most significant in affecting the length of stay 
in shelter. 
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TABLE VI 


ll 

/ 


/ 


Sensitivity Analysts of Time Required In Shelter 


Parameter 

1 

Mean Value 

*1 

Gradient 

Sensitivity Index 
(Equation 21) 

Sf 

Z 

1.260 

- 144.000 

-3.860 

F 

0.914 

- 144.000 

-2.220 

‘^3 

144.000 

0.373 

1.140 

B 

0.025 

-1700.000 

-0.918 

RP 

79.600 

0.331 

0.560 


The Independent variables (which affect the length of stay in shelter) are 
ranked on the sensitivity Index In decreasing Importance as follows; 


1. Field decay exponent, Z “3.9 

2. ERD recovery fraction, F • 2.2 

3. Maximum ERD, » 1.1 

4. ERD recovery rate, B ■0.9 


5. Reference intensity/protection factor ratio, RP "0.6 

Of equal interest to the parameter sensitivity analysis is the actual 

timing of Maximum ERD and the number of days' stay in the primary shelter before 

recovery to an allowable dose Is accomplished. For a shelter PF of 73, the 

Maximum ERD of 144r occurred on the 7th day after the attack. On the average, 

3/ 

47 days elapsed before the ERD recovered to the allowable dose level of 80r.— 
Figure 9 shows the frequency of occurrence plotted against the duration in days 
before reaching an ERD of 80r. It can be seen that there are cases in which 
stay time would exceed 100 days if the ERD of shelter occupants is to return 
to 80r. 

K, Evaluation of Significance of Sensitivity Indices in Shelter Systems Analysis 
1. General 

Hie objective of shelter systems analysis is to predict as accurately as 
possible the adequacy of existing or potential protective measures against 

3/ 

— If instead, an allowable level of 50r (Subgroup 19) were chosen, then the average 
length of stay would increase by 33 days (from 47 to 80). 
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Days for Return oi.' ERD to 80r (N^) 

Fig. 9. Frequency of Occurrence of Stay Times vs Days for Return of ERD 
(for Cases in Which Max ERD < 200) 










fallout. In general, the better the parameters that determine the adequacy are 
knoim, the better the predictions will be. 

The sensitivity indices derived in the above analysis give an indication of 
the effects of .^ot knowing the various input parameters exactly; the analysis 
does not yet take Into account how accurately these parameters might be known. 

In the following sections, an attempt will be made to Identify the approxi¬ 
mate magnitude of uncertainty to be expected in the eight parameters considered. 

2. Reference Intensity and Time of Ar-ival 

Of the eight parameters that were considered in the sensitivity analysis, 
reference intensity and time of arrival are largely dependent upon the choice 
of strategy in the attack; consequently, these parameters cannot be known 
accurately. Desired ground zeros and the occurrence of air or surface burst are 
the choice of the attacker. 

At any given point, the reference intensity and time of arrival are also 
affected by wind speed and direction. 

The RISK type programs are attempts to represent the combined effects of 
strategic decisions and wind conditions. In the choice of input variables for 
the sensitivity analysis (section II.A), it was seen that reference intensities 
between 300 and 3900r/hr and times of arrival between 1 and 7 hours seem to 
represent the range of intensities with which fallout-only shelter systems are 
primarily concerned. Variations from the means of these input parameters of 
+ 75% would seem to be a representative estimate of the probeble range of values 
for reference intensity and time of arrival in an analysis of a shelter system. 

3. Buildup Factor 

The values of the buildup multiplying factor used in the sensitivity 
analysis ranged from 1.13 to 6.15. These extreme values vary from their mean 
by about 70% and this is proposed as a range of values that could be expected 
in systems analysis. 

4. Protection Factors 

The National Fallout Shelter Survey (NFSS) produced data on the present 
number and the protection factors of shelter spaces. However, there are various 
reasons why these data alone may not represent the actual shelter posture of the 
population. 
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First, based on a limited sample of shelter [Reference 5], it appears that 
procedures used in the Phase 1 shelter survey generally underestimated the 
actual protection factor by a considerable amount. The current updating of 
NFSS data is graduaHy improving the accuracy of information on shelter status, 
but a complete rcevaluation of all NFSS shelters is not contemplated by OCD. 

At the present time the average potential protection factor of the United States 
population may be considerably more than presently represented in Phase 1 or 
Phase 2 data. However, ever, if the protecttur. factors wvre accurately deter¬ 
mined, there still would rem.-'in uncertainties as to shelter utilization that 
would influence the results of shelter analysis. The question of shelter 
utilization was presented in Volume II where it was considered in terms of 
actual SMSA's. (See Volume II, A Sensitivity Analysis of Selected Parameters 
Based on 8 SMSA*s.) 

Consequently, random variations of 80% about the mean can be considered 
reasonable in estimating the significance of the sensitivity index for protec¬ 
tion factor. 

5. Time in Fallout 

Estimates of length of time outside in fallout are arbitrary. By assuming 
adequate warning, no one would need to be outside. However, there may also be 
cases where large numbers of the population are in fallout for periods ex¬ 
ceeding the 0.6 hours arbitrarily used in the sensitivity analysis. 


Although the sensitivity index for time in fallout is small, the importance 
of this parameter in systems analysis must be determined in conjunction with 
actual locations and categories of shelters. The combined effects of these arc 
analyzed in Volume II. 

6. Decay Rate - --- 

The decay rate may vary with the construction of the weapons and to a lesser 
degree with the nature of the surface materials at ground zero. Variations 
probably will not be more than plus or minus 25 percent of the 1.2 rate of 
decay, which is based on empirical data. 

7- Biological Factors 

The fraction of the dose that is reparable (F) and the rate of recovery 
(B) are important parts of the Blair Formula for calculating ERD. However 
according to Reference [6], a firm experimental basis appears to be lacking 
for the values presently used. The author in the same reference also questions 








whether the constants are Independent of time. 

In this study, the average value of F, the recovery fraction, is taken as 
0.90; variations of 10 percent would seem likely. The average value of B, the 
recovery rate, is 0.025 and probably could vary by 70 percent. 

8. Comparison 

We then may compare the sensitivity Indices with the rough estimates of the 
accuracy to be expected in input variables. These are sunnarized in Table VII. 
An approximation of the relative importance of each parameter in contributing to 
total uncertainty in shelter analysis is obtained by multiplying the range of 
uncertainty by the sensitivity l.-^dices. 

TABLE VII 

Comparison of Sensitivity Indices and Uncertainty in Input Parameters 


Variable 

S (ERD 
Model) 

S (Total 

Dose Model) 

Uncertainty 

Approximate 

Relative 

Importance 

Field Decay Exponent 

-4.22 

-4.36 

25% 

1.09 

Refererce Intensity 

+1.02 

+1.00 

75% 

0,75 

Protection Factor 

-0.96 

-0.99 

80% 

0.80 

Time of Arrival 

-0.50 

-0.50 

75% 

0.35 

Buildup Duration Factor 

-0.18 

-0.25 

70% 

0.18 

ERD Recovery Rate 

-0.14 


70% 

0.10 

Time Outside in Fallout 

+0,05 

+0.02 

100% 

0.05 

ERD Recovery Fraction 

-0.30 

--- 

10% 

0.03 



CONaUSIONS AND RECOMMENDATIONS 


A. Conclusions 


Considering both Che sensitivity of each parameter and the range of uncertainty 
Chat each is likely to have, it is possible to draw some conclusions as to their 
effect on Che results of systems analysis. 

In analyses for estimates of Maximum ERD, Che following have t!ie most signi¬ 
ficance: Field decay exponent (Z), reference intensity (I^)* and protection factor 
(PF). Having less significance are time of arrival (T^) , recovery fraction (F), 
ERD recovery rate (B), and buildup duration factor (E). 











Time in fallout (T 2 ) has a low sensitivity index. However, extending the range 
of values that were assumed for it beyond 0.6 hours may have a considerable effect 
on the magnitude of this number. The importance of time in fallout was better shown 
in Volume II, where it was analyzed as a factor in shelter utilization. 

In analysis of required time in shelter, field decay exponent (Z) and ERO re¬ 
covery rate (B) are the most significant. 

of the above parameters, reference Intensity and time of arrival are almost 
entirely dependent upon the strategic choice of the attacker and wind conditions. 
Consequently, better estimates of their effects in systems analysis is dependent 
upon military intelligence more than upon additional research and analysis. 

The importance of the field decay exponent in determining both fatalities and 
shelter stay mkes it important to have continuing analysis of the validity of the 
1.2 value. There may be appreciable gains in accuracy possible since Che exponent 
is primarily dependent upon Che physics of the radioactive decay process. 

Increased knowledge of the shelter protection factors available to the population 
is important for Che making of valid systems analyses. The effect of underestimates 
of protection factors can be taken into account by using a multiplying factor. 

However, it is Important to define the factor more closely than is possible now. 

This would not necessarily require large expendit’.res for obtaining new data or for 
corrjccing data Chat now exist. Statistical analysis based on Che parts of the 
survey that have been updated or examined in detail could provide a much better 
estimate than evidently now exists of the overall shelter that would be available to 
the population. 

Identification of the biological recovery process is particularly important In 
estimating time required in shelter. Possible approaches to reducing tne range of 
uncertainty of this input are given in Reference [6]. 

Finally, efforts to define more closely the limits of any input variables should 
be viewed in the light of their relative effect on output data. For example, attempts 
to define protection factors accurately may make a change in the relative importance 
of that parameter that is small compared to the relative importance of reference 
intensity, which cannot be accurately predicted by any research means. 

The sensitivity index of the buildup duration factor (E) is based on "moderate" 
conditions where the time of arrival approximates 3.3 hours and (E) approximates 2.6. 

In observing the value of the sensitivity index S^, it was found that the index is 
quite dependent on changes in (E). Some further analysis of the effect of the build¬ 
up duration factor is discussed in Appendix A. 
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Ihe analysis based upon return of ERD to 80r (the threshold of appearance of 
clinically detectable symptoms of radiation sickness), showed that shelter stay times 
of much more than two weeks could be necessary. It should be noted that this analysis 
included only cases in which the Maximum ERD was 200r or less. Without these restric¬ 
tions, even longer required stay times would result. This is very significant in 
planning for the stocking of shelters, recovery planning, and for the establishment 
of policies for postattack management of supplies. 

In summary the conclusions from the sensitivity analysis are: 

(1) Casualty calculations are quite sensitive to errors in the field decay 
exponent. Ihe sensitivity of dose or casualties to the field decay cate 
remains at approximately -A.O over the range examined. (The decay rate 
was varied from Z «• 1.1 to Z = 1.4.) 

,(2) Sensitivity to variations in fallout reference intensity and protection 

factor are high.over the whole range of parameter values. Hence, precise 
knowledge of the fallout shelter posture and the fallout reference inten¬ 
sity is much more essential to accurate vulnerability analysis than the 
remaining parameters in most cases of interest. 

(3) Sensitivity to time of arrival of fallout can be quite high, in some 
radiological environments. 

(4) Sensitivity of casualty computations to the remaining parameters is low 
in most cases of interest. 

(5) Expansion of the sensitivity analysis to include parameters other than 
fallout, which define the total casualties from a given attack on the 
United States, is necessary before further conclusions concerning a 
national vulnerability analysis can be drawn. 

B. Recommendations 

1. Vulnerability analyses should employ protection factors computed by best 
available methods, and research and/or surveys to improve protection factor data 
should be encouraged. 

2. Because or the sensitivity of systems analysis results to the field decay 

-1 2 

exponent, continuing analysis of the validity of the t * decay law should be made. 

3. Additional study is required to establish the sensitivity of fatalities, 
casualties, and dose to duration of shelter stay. 
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A. Sensitivity analysis should be extended to include the parameters defining 
the effects of blast and alternative measures of effectiveness (casualties by type, 
lose distribution of survivors, etc.)* 


5. Sensitivity analysis should be applied to identification of the important 
cost/effectiveness parameters used in the budget allocation model (see Volume I, 

\ Cost/Effectiveness Computer Procedure for Optimum Allocation of Fallout Shelter 


System Funds Under Uniform or Variable Risk Assumptions' 
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Appendix A 


The Intensity Buildup of a Radioactive Fallout Field 


In calculating the radiatior dose received at a point in the fallout field, it 
is necessary to determine the rate at which the fallout particles are deposited from 
the'cloud. The easiest method of representing this buildup is to plot the fraction 
of total fallout mass deposited i.s a function of tire. If is the time the first 


parp|.cles arrive, and is the lime the last particles are deposited, then the 

duration of the buildup period, - T , is approximated by a relation of the form 

C . A 

■ '^A ’ ^A • " 

where E is the buildup duration factor referred to in this volume. 

I 

Figure A«1 presents the buildup function [Reference A-l[ which was obtained 
from weapons tests. In addition^ two approximations to this function arc plotted 

I 

in F^igure A-1 for = 1 and E *=,2.5. The analytic form for these approximate 
curves is given by Equations A-2:and A-3. 


j<'> ■ i‘e • (' ■ e) 

F(t)| . i - cos 


<A-2) 


(A-3) 


Equation A-3 gives a somewhat better fit of the curve fi-om Reference '^A-1', 

[ ’ * . 
hence this buildup function was used in the computer model that calculated Maximum 

r 

ERD. However, the total, dose expressions cannot be integrated in closed form if 
Equation (A-3) is used, hence Equation (A-2) was used as the buildup function for the 

analytical total dose model. As might be expected from Figure A-1, little difference 

' 

was found between sen.sitivity indices from the two models. 

j REFERENCE 

j . 

[Afl] Miller, C. F. Fallout and' Radiological Countermeasures , SRI Project IM-4021, 
Vol. I, January 1963. 
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Appendix B 

_ The Mainline Model and Computer Program 

I, THE OPERATIONAL TRANSFORM OR MAINLINE MODEL 


Attack Environment Formulation 


1. Introduction 

The first step in formulating the Mainline Program is proper characteriza¬ 
tion of the attack environment. The interrelation between the Initial intensity, 
field buildup, the field decay, and the field attenuation by shelter, are the 
four functions that predict the dose rate. To determine their relative im¬ 
portance, it becomes mandatory that tnis interrelation is accurately formulated. 

The underlying assumptions throughout this section are: 


(a) The Reference Intensity of the radioactive field (in roentgens per 
hour) will conform to the measurements given in the NREC RISK-type 
data (i.e., an intensity (I^) at H + 1 hours) as well as the distribu¬ 
tion of intensity and arrival of fallout obtained also from NREC RISK- 
type data. I 


(b) The measure of protection obtained by the emploment of shelters will 
be the Protection Factor (PF) obtained from ’'Th«i| 

Phases 1 and 2. 


I 

Fallout Shelter Surveys 


2. Decay of Field Intensity 


Since a radioactive particle decays with time, our problem is to determine 
the appropriate mathematical representation, when applied to residual fields 
used in Civil Defense Systems Analysis. The widely used function shown in 
Figure B-1 and Equation (B-1) will be used. 








where 


f)(t) ■ Field Decay Function, in r/hr. 
■ Reference Intensity, in r/hr. 
t ■ Time in hours. 

Z *■ Empirical Constant. 



Fig. B“l. Decay of Radioactive Field 

The two important constants introduced with this function are: the 
empirically derived exponent (Z), and the reference intensity (1^). The refer¬ 
ence intensity is one of the two interface parameters connecting the fallout 
shelter analysis model with the 'S/ar games" which establish the probable attack 
environment. 

3. Buildup of Field Intensity 
a. Introduction 

The problem now is to characterize analytically the field intensity 
during the buildup phase. The fact that there is no existing data within 
our files on which the usual methods of curve fitting can be applied, forced 
us to use a more general form of analysis. Our analysis uses what little 
significant information there is and employs as much intuition as seems 
necessary, useful, or feasible. Although there are two relevant functions 



(B-2) 


and 




(B-3) 
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they are incompatible over a wide range of times of arrival. Tlie relation 
time of peak intensity equals twice the time of arrival--is taken as being 
more nearly valid. In the field buildup, two distinct phenomena are 
present. The first is the arrival (at ground level) of the radioactive 
particles. Hie second is the independent phenomena of the radioactive 
decay of these particles. Thus, characterizations of these two phenomena 
will be independent, and the field Intensity represented by the product. 

The hypothesis implies a particle buildup function of the following 
general form given in Figure B-2, which, when integrated, yields the time 
dependent arrival function given in Figure B-3. 



Time t -► 

Fig. B-2. Particle Arrival Rate 



Fig. B-3. Accumulation of Particles 
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The product of the functions represented in Figure B-1 and Figure B-3 
would give the field intensity shoim in Figure B'4. 



Fig. B-A. Field Intensity vs. Time 

b. Formulation 

From several general forms a displaced Cosine function was chosen 
for simplicity and compatibility with the function (t •• 2T.) to represent 
the fallout particles accumulated (Reference Equation B-A). 

F(t) - J - cos (— 

where 

T, < t < T„ 

A — — C 

and 

1 ■ ■'c • ■'a- 
-z 

The field decay function (i.e., t ) times the buildup function times the 

reference intensity yields the final representation of the field Intensity 

over time. For time between T. and T : 

A C 

“ cos(- 7 ^ ")] (B-5) 
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Thus, the characterization of the field intensity is in two parts: the 
first during accunulation--fron (time of arrival) to (time of 

cessation). The second, after all fallout has occurred (Reference Equation 
B-1), describes the field. 

In order to comply with T = 2T., the value of n must be fixed. The 

p A 

following are values of for three values of Z: 

Z - 1.0 T] - 

Z - 1.2 r. - 1.175T. 


Z - 1.4 


1.20T. 


Attenuation of Field Intensity 

a. The Problem 

The population or an individual may move through a series of shelter 
conditions over time. The duration of each shelter phase, although 
nominally under operational control, has well-defined bounds based upon 
the attack environment. The initiation of the recovery phase depends upon 
minimizing the stay in shelter. Thus, the influence of shelters when 
modeled must take into account the operational requirement to minimize 
shelter stay time. 

b. Formulation 

Fallout shelters are classified by a protection factor (PF). Tlic 
protection factor function operates on the unprotected radiation field 
to yield the attenuated field that exists within that shelter. 

There are four major operational phases (Reference Figure B-5). 

These are: 

(1) An initial unsheltered condition which corresponds with a person 
moving through a fallout field before entering shelter. This 
phase will be called "time outside."—^ It is related to the 
efficiency of warning systems, drills and training, and shelter 
assignment plans. 

(2) Phase 2 is the primary shelter phase. In most current models it 
is the only one considered. Operational plans usually specify 


—^ The symbol used in this report for initial time outside in a fallout field is T 2 . 


B-5 





Field InCensiCy (r/hr) 


( 3 ) 


Phase 3 is operationally assoc 
begin recovery. The objective 


a two-week stay, and they are Stocked on that basis. The control¬ 
lable variable is the length o^ stay in the phase. Our model 
considers the length of stay dependent upon an objective of the 
third phase and is subject to explicit calculation. 


ated with emerging from shelter to 
in this period is to obtain produc¬ 


tive effort (usually at a PF or 1) without any physiological 
radiation symptoms. To do this introduces a decision parameter 
(an allowable ERD), an operational parameter, and the equivalent 
PF required. The allowable ERD lies between the limits of 
clinical effects (25r) and symptomatic effects (lOOr). 


(4) Tne fourth phase is the return 


to normal. It begins when the 


i|idividual's recovery rate equals the dose received in the un¬ 
sheltered field. It is characierized by the lack of control 
variables, decisions required, etc. In other words, the effects 
of normal biological recovery processes exceed the effects of the 
radiation fields. 

The problem domain investigated in this volume extends only to P’hase 3, 
time (Tj). (Practical considerations precluded extensive investigation of 
Phases 3 and 4 under the present contract.) 



Fig. B-5. Attenuation of Field Intensity 
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In Fig. B-5, the solid line represents the sheltered point intensity. 
In Fig. B-5, the dotted line represents the field Intensity if: 


where 


PFj ■ PF^ ■ PFj - 1. 


Phase 

1 

PFj^(t) - 1, 0 < t < Tj, 

Phase 

2 

PF2<t) - PF, T 2 < t < T^ 

Phase 

3 

PF^(t) » Equivalent PF(t), T^ < t < T^ 

Phase 

4 

PF^(t) - 1, T^ < t. . 


5. Summary Attack Environment Model ' ' ... -.. 

The attack environment radioactive field intensity at any time (t) is the 
product of: the Reference Intensity (1^)» the particle buildup function (F(t)); 
the field decay function ('(t)); and the protection factor function (PF(t)). 


«(t) - R . F(t) 


9(t) 

PF(t) 


(B-6) 


Or by substitution: 




when T, < t < T, + -|T 

A ^ A A 


t(t) -0 when t < (B-7) 

I. 

t'Ct) - t when t > + r)T^ 


where: 

ti(t) ■ Field Intensity in r/hr at time t after the originating 
explosion event. 

I^ = The reference intensity in r/hr, one hour after the explosion 
event. 

t ■ Time in hours measured from H. 

T. ■ Time of arrival measured from H hour. 

A 

PF, ■ Protection factor at time t. 
k 

Z ■ empirically derived field decay constant. 

n ■ T_ - T, ■ fallout interval. 

* C A 
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Biological Response 


1. Introduction 

Ihe state, efficiency, and outcome of a shelter system analysis is measured 
by the biological Impact, the usual criteria being the number of fatalities or 
lives saved. Fatalities are calculated by using a probability of response to 
the total radiation dose received. 

At the present state of understanding, the equivalent residual dose model 
is the best available formulation of biological response to dose received. 

2. Equivalent Residual Dose 

To represent analytically the absorbed gamma radiation field in terms of 
a time dependent dose "state" requires taking into account the observed phenomena 
of biological recovery. The ERD model relates the field intensity and biological 
recovery rate to give a single "dose" value, approximating the physiological 
state of a human. 

The ERD at any point In time consists of two parts. The nonrecoverable 
portion is the accumulative amount of the nonrecoverable dose received over 
time. The recoverable portion has two additive parts; one is a fraction of the 
field contribution within the unit time considered, and the second Is the 
recoverable dose carried over from the previous period. The basic statement 

of ERD (recoverable dose (D ) in time increment, n) is the sum of that not yet 

^ th 

received from previous periods and the new contribution during the n time 

period. 

D^ - /" (B-8) 

n-1 

where ■ __^____ 

B is the percent per unit dose recovered in period n, and F is the 
recoverable fraction. 

Nonrecoverable dose (d ) in time increment (n) is the accumulative sum of all 

n 

nonreparable doses: 

,dn - d^_j + (1-F) /" <5(t)dt. (B-9) 

n-1 

The practical problem of describing the function of the ERD(n) in a closed, 
well-behaved form requires that it be approximated by finite differences. 

Finite difference problems, which require recursive solutions, lead to computer 
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implementation. A solution is to use the basic difference equations with n 

in one day intervals and with d being computed daily, then added for the 

^ th 

equivalent residual dose on the n*^ day. 
where 

D ■ Recoverable dose in time interval n-1 to n, 
n ’ 

d ■ Nonrecoverable dose in time interval n-1 to n, 

■ Reference intensity, 

F <■ Fraction of dose which does reparable damage, 

(1-F) ■ Fraction of dose which does nonreparable damage, 

B ■ Dally fraction of dose repaired, 

-■ The time of arrival of the first measurable radioactive fallout, 

Z ■ Decay constant, 

PF ■ Protection factor. 

3. Probability 

Ihe problem is now to relate the equivalent residual dose which an individ 
ual accumulates to the probable biological state. The usual specifications for 
radiological induced states are fatalities and casulaties, with fatality being 
a type of casualty. Because this causes some confusion, we will introduce the 
term noneffective. For a given Maximum ERD, we will then have three probabili¬ 
ties: 

(1) Death. 

(2) Noneffective--The condition of being alive but physically unable to 
perform a productive task. 

(3) Normal and Marginal Effective—Individuals with acceptably low Overt 
physiological symptoms (below lOOr). 

The following linear functions are used to approximate the empirically derived 
probability of casualty vs. ERD curves (reference Figure B-6). 

a. Probability of Fatality 

Pj - 0 when ERD < 200r 

. vhen 200r < ERD < 700r 

Pj « 1. when 700r < ERD. 
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(B-10) 









g lOOr 200r 300r 400r 500r 600r 700r 

4 * 


Equivalent Residual Dose (ERD) 
Fig. B-6. Dose Response 



Probability of Noneffectiveness 



when ERD < lOOr 
when lOOr < ERD < 200r 
when 200r < ERD < 300r 
when 300r < ERD < 700r 
when 700r < ERD. 



(B-11) 


(B-12) 


C. Mainline Program 
1. Introduction 

Ihe models, previously described, are combined Into a single msdel—a 
computer program. The finite difference mf^thod, particularly adaptable to 
computer prograrming, is necessary because of several functional discontinuities 
which preclude integrating the individual models simultaneously. In translating 
algebraic formulas Into FORTRAN coding, a change in notation Is required. To 
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assist in identifying corresponding Ct.rn.s, a Clossary of S>T.brls osc-J in the 
progransning, and the corresponding syrbols found in ratherarical for;-iulatiori, 
are listed in Table B-I. 


TABLE B-I 
Progran Syr^bols 


Parameters 

Model Symbols 

Pingraming Symbols 

Reference Intensity 


K 

Time of Arrival 

. ^A 

I/. 

Time Outside 

^2 

•f ^ 

. >. 

Allowable Dose 

^4 

!>■ 

Field Decay Constant 

Z 

Z 

Recovery Fraction 

F 


Recovery Rate 

« 

9 

Buildup Duration 


r.'A 

Protection Factor 

PF 

PF 

Protection Factor Categories 


K 

Maximum No. PF Categories 


KM.V; 

Time (the variable of 
integration) 

Hours 

Days 

V 

w 

n 

T 

Total (ERD) Dose 

EBD (n) 

D 

Recoverable O-'se 

D 

n 

D1 

Konrecoverdble Dose 

d 

n 

D2 

Day Maxirrjm Occurred 

• T 

N1 

Field Intcrnity 

:(t) 

RDN or FLr 

Day Allcrk'able Dose Occurred 

^2 

N2 

Probabiiity of Fatality 


PI 

Pribability of Noneffccrlve 

^2 

P2 


In the input sections of the progr.in, if a programing symbol given in 
Table 3-1 is folioved by L, H, or S, the ruirbvr representt-.l by the symbol in 
either the lowest, higi^est, oi In.crevnentai vala taken by the variable repne- 
aented by the pre'.aJing letters (i-t., TAL rs.'ana the lowest value the time of 
arrival will rake). , , 


R-ll 








tI *re '*•? fl** ii*sr*T! by of|^.i*»ticn' cf 

tbf •tatr-rw-rt lire rii«*-et vith tb« . Specifically, each of the 

codlof lire rariefi art four digi ' - • lr«t digit representi the najor routine 

being calculated (tee Table B-II); i.e aecond and third digits correspond to the 
flow dlagran box nunber (xx); the fourth or unlts-posltlon digit Is "open" to 
allow for sequencing corrections, additions, and expansions (0). 

TABLE B-n 

. .. _ Major Program Routines ______ _ __ J 

Code 

Number Series 


Input 

IxxO 

- - Iterations & Loops- 

--2xx0 - 

Calculation of Buildup 


Duration in hours 

3xx0 

Initial Conditions 

4xx0 

Calculate 1st Day Doses 

5xx0 

Day by Day Doses 

6xx0 

Outcome Probabilities 

7xx0 

Formats & "Housekeeping" 


Steps 

8xx0 

Explanation of the Flow Diagrams ^ 


The first section of the program (block numbers 1-36) Initializes the 
input data and sets up the control loops for the variables. Block numbers 36 
through 40 select the appropriate duration multipliers to calculate the dura¬ 
tion of Particle Buildup (T3), 

T3 - (1.133 + ETA)*TA when Z =1. 

---*^(17t75 •h-ETA)*TA-whenrZ^irivZ-^-—^ (8^13) 

• (1.200 + ETA)*TA when 2 * 1.4 

the range of Buildup Duration Factor (ETA) Is between zero (corresponding to 
the displaced sine function) and five (corresponding • .'th 5TA ^ - TA function). 

The program sets the initial conditions In block 47. Clock Time (T) begins 
at time of arrival (TA) and the Day Counter (N) begins at one, the time is 
Incremented by 3 minutes or 0.05 hrs. 
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Field ir.tcrsity cir -jnit jf tirif: ii calculated ir, blocVs i-, i5, 

uiinR the product of Reference irtensity (rA<r), Particle Arrival functlcr, arc! 
Decay function. 

Llnlt on T1 Is: 

RDN - I jl - cos(^^)sJt”^ 0 < < r. (B-14) 

Ihe test, block 42, determines if there is any movement outside in fallout-'or 
if T2 is positive. 


Block 52 is similar to 47 except the field is attenuated by the protection 
factor. Block 53 tests to determine if buildup of particles is continuing. 
Block 54 increments time out of buildup phase so that field intensity is 
governed only by the decay function corresponding to the algebraic formula. 

(5 - . (B-15) 


Block 56 tests to determine if the 1st day is completed. Block 55 increments 
time by 12 minutes during buildup. Blocks 49, 50, 51 are identical to 44, 45, 
46 but are not linked because of a change in the time intervals. 


FLD in block 59 is equivalent to RDN but is used to account for dose when 
th>. time step is by days. 


-(t) 


.24.HA 


,1-2125+14 

r(t) • It when Z i* 1 (B-16) 

CM. iT*:.LuC2^tiAii'ii 
'' ' Z-1 


Block 60 does not have a time term since the integration in block 59 was over 
a finite time. “ 


Blocks 64, 65, 66, 67 and 68 sets up transition from the first day, to 
day-b>-day caicui^tions. At blocks o4-C3 tirau ERD division into recoverabl 
and nonrecoverable fractions is made. Block 66 steps time by 12 hours. Block 67 
sets maximum dose register. Block 69 solves equation for Field Intensity summed 
over 24 hours but corrected according to PF, 


__ R*t"^*( 24 hrs.) 

-PFOO-^ 


(B-17) 
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T2 ti t TtiovtTdoit, Biock 71 «cCu?sblAtc* ttie ncfiiecy/erai !< 

fioif. Si«KU 72 teit» for •arc liearh. Block 74 test* the «!«y'* <lo»e vt. 

E310. If there It a new Maxirun EKC, it !• let In block 75, and the day and £1.-% 
arc increnented In blocks 76 and 77. If t^'e peak ERD (D3) has been passed 
(D < D3), then the next point of interest is whether or not the allowable level 
has been reached. This test is block 78. If the dose is above the allowable, 
block 79 testa to see if the allowable is reachable--!.e.,nonrecoverable dose 
is less than the allowable. If the nonrecoverable dose exceeds the allowable, 
then the program is in a loop as the length of stay in shelter would approach 
infinity—therefore, the error message. 

Blocks 80 through 92 calculate the probability of fatality and the probabili 
ty of noneffectiveness. 

The preceding explanation should give some understanding of the computer 
program and how it represents the mathematical formulation. This program, then, 
generates the possible data combinations that are analyzed in this volume. 









/ ^ 
i 
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Code Kuinber Series IXXO 
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) 

SET 


Total ERD 

D=0. 

Recoveralle dose 

D1=0. 

Non-rec. dose 

D2=0. 

Max, ERD 

D3=0. 

Day under 


calculation 

N=1 

Day max. ERD 


occurred 

Nl-l 

Day ERD return 


to allowable 

N2 = l 

Time in hrs. 

T=TA 


@r 


Increment Time 
by 3 min, 

3;_»_T + .. 05 .. 


mm 

■RSS 


‘ , 13 . 

1 _ 

_ 


©f 


RDN 

= 1. - C0SF(RDN)j 


RDN 

_ R*R0N*T* 

2.0 


= D + RDN*.05 



Code Number Series 4XX0 
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Code Number Series 6XX0 
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PKO&KAM MAINLINE 

DIMENSION PFNdOI ___ __ _ .. 

N2-N2 

KO0£*0 - _ ___ 

READ INPUT TAPE StlOUb.RL»RH,RS 

_ READ INPUT TAPE S.lOQQ.TAL.TAHtTAS __ 

READ INPUT TAPE 5,1000,T2LtT2H,T2S 

REAO^INPUT TAPE 5, lOOOj Q4L • D4HiP_45_ __ 

READ INPUT TAPE 5t1000«ZL.2H.2S 

READJNPUT TAPE 5. lOOO.FL. PH.F S __ __ ^ 

READ INPUT TAPE St1000«BL,8H,8S 

_ READ INPUT TAPE S» 1 QUO. ETAL .ETAH. ETAS _^__ 

READ INPUT TAPE 5.lOQl.PFNIl11PEN I 2)tPFNI 3 11PEN(AI,PEN(S11PFN(6), 

lPENm,PFNI8)iPFNiq)__ _ 

READ INPUT TAPE SdOOZtKMAX 

1000 EJDRMATI 3E15d)__;_ 

1001 EdRMAT(9F6.0) 

1002 FORMAT! 110) _;_^___ 

MKITE OUTPUT TAPE 6»89S3 

8933 FORMAT! 114H_R_TA_12__DA,_ I __F Q 

1 ETA K Pi P2 D3 N1 0 N 02 I 

2110 _RfRlrJiS __:_:... 

2120 i[F(K-RH)2130.2130.8820 

2130 H«R»RS __ 

21A0 TA»TAL-TAS 

2150 lF(TA-TAH>2160t2160.212 0 _____,_ 

2160 TA=TA*tAS 

2170 _T2*T2L-T2S_ __ . ____ 

2180 IF IT2-T2HJ2190,2l90t215b 

2190 T2 = T2»T2S __^__ 

2200 DAsOAL-OAS 

2210 IF (OA-OAH>2220.2220.2160__ 

2220 DA»bA+OAS 

2230 2-ZL-ZS ___ 

22a 6 lF(Z-ZHi2256.2250.2210 

2250 __ 

F=FL-FS 

2 2 7 0__ 1 F(F-FH) 2 28 0.22 80.22 A0_____ 

2280 F=F+eS 

' _B»BL-BS_j___ 

2300 1F(B-BH)2310.2310.2270 

2310 B=B»BS __ ' __ 

ETA=ETAL-ETAS 

2330 IFIETA~ETAH)23A 0.23A0.2300 __ . __ 

23A0 ETAxETA+ETAS 

2350 DO 9000 Kxj.KMAX_ ^ __ 

3360 IF I 2-1.113370.3380,3380 

3370 T3x( 1.1333*ETA)»TA __ 

GO TO AOOO 

_3380 1F<Z-1.13)3390.3A00. 3A00 ___ 

3390 f3=l1.175*ETA)»TA 

__ GO To_Aopo_:_ _ 

3AbO "T3 = ( 1.2b*EtA)*TA 

AOOO CONTINUE _ 

AAIO D»0. 

pl-0. _ __,_ 
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6620 

T»TA 

IF(T-(T2+TA>16630,5680,5680 


6630 

T-T*0.05 


6660 

RDN-(T-TA)*3.i6l6/T3 


6650 

r6n*i.o-cosf(ron> 


_ 6660 

R0N»(R»RDN>#(T»*ZZ>/2.0 


6670 

0«04'(R0N*0.05> 



GO TO 6620 


5680 

T»t*0.1 


5690 

R0N»|T-TA)*3.1616/T3 


5500 

R0N«1.0-COSF(R0N> 



R0N*(R#H0N>»(T«»ZZ>/2.0 


5520 

0*0+(R0N*0.2/PFN(K> > 


5530 

IF(T-(TA^T3>>5550.5550.5560 


5550 

T«T^0.2 

GO TO 5690 


5560 

T«T»0.1 

IF(T-(TA^26.0>>5570.5610.5610 


5570 

IF(2-1.0>5590,5580,5590 


5380 

ARG=26.0/T 



FL0>L0GF(ARG> 

GO TO 5600 


5590 

ZZZ*1.0-Z 

FL0*l(T*»2Z2>-((TA*26.0>**ZZZ>)/(Z-1.0> 


5600 

D«0+(FLO»R/PFN(K>> 


5610 

IF(0-700.>6000.9000,9000 


6000 

C(3NTINUE 


6660 

D1»F*D 


6650 

02»(1.0-F>*0 



66 70 
6660 


6730 

6760 


6750 

6760 


6780 

6790 


7710 

7000 


7800 

T9 1 Q 


7820 

7830 


03«0 

N»N*1 


6690 

6700 


6710 

S720 


02»02*< l.pl-FUFLO 

■D«fll*P2 


FL0*IR*26 ,!0»«CT**ZZI/PFN(K» 
01«(l.-BIP01*(F*FLD) _ 


IF(0-700.16760,9000,9000 
IF(0-0316780.6750,6750 
03-0 


N1 


6770 T-T+26,0 


GO TO 6680 _ 

IF I 0-06>7000,6790,6790 
IF(02-06>6770,7710,7710 
KOOE«1 


KOOE«1 
CONTINUE 


N2«N 

IF(03-1Q0.>9000,7820.7620 


IF(03-200.I 7880.7830.7830 
IF(03-300,>7890,7860,7860 










7B40 IF(U3-700.)7900,9000,9000 

7880 P2 = <03-100. )/?00. _ ___ _ .. 

7910 Pl=0. 

CO TO 8000 _ . _ ___ 

7890 P2=(tJ.O*D3»-100.J/1000. 

_ CO TO 7920 __ 

7900 P2=l70n.-n3>/500. 

7920 Pl*(03-200,1/500. . _ .. __ . 

8000 CONTINUE 

8950 V(RITE OUTPUT TAPE 6, B934j KOOE » R, T A, T 2,04, Z , F , B, ET A, PFN ( K ), 
iPl,P2,D3.Nl,0,N2,02 

8934 FORMAT!I4.4X.F6.0.FH.l.Fe.i.F6.0.F6.2,2X,F5.3,2X.F5.3.2X.F5.3. 
1F5.0,F7.3,F7.3,F7.1,I5,F7.1,15,F9.2> 

9000 CONTINUE __ __ ___ _ . 

KODE»d 

__ _G0 IP_2310 ___• 

8820 CONTINUE 

.END SENSPROG .■ .^ .-.....- ....-.. 


LNO 

PLNl^ 


300. 

3000, 

1200. 

I.O 

6.0 

3^0 

d.u 

.5 

.3 

20. 

70, _ 

30. 


1. 1,3 
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